Extraction-free whole transcriptome gene expression analysis of FFPE sections and histology-directed subareas of tissue.

We describe the use of a ligation-based targeted whole transcriptome expression profiling assay, TempO-Seq, to profile formalin-fixed paraffin-embedded (FFPE) tissue, including H&E stained FFPE tissue, by directly lysing tissue scraped from slides without extracting RNA or converting the RNA to cDNA. The correlation of measured gene expression changes in unfixed and fixed samples using blocks prepared from a pellet of a single cell type was R2 = 0.97, demonstrating that no significant artifacts were introduced by fixation. Fixed and fresh samples prepared in an equivalent manner produced comparable sequencing depth results (+/- 20%), with similar %CV (11.5 and 12.7%, respectively), indicating no significant loss of measurable RNA due to fixation. The sensitivity of the TempO-Seq assay was the same whether the tissue section was fixed or not. The assay performance was equivalent for human, mouse, or rat whole transcriptome. The results from 10 mm2 and 2 mm2 areas of tissue obtained from 5 µm thick sections were equivalent, thus demonstrating high sensitivity and ability to profile focal areas of histology within a section. Replicate reproducibility of separate areas of tissue ranged from R2 = 0.83 (lung) to 0.96 (liver) depending on the tissue type, with an average correlation of R2 = 0.90 across nine tissue types. The average %CVs were 16.8% for genes expressed at greater than 200 counts, and 20.3% for genes greater than 50 counts. Tissue specific differences in gene expression were identified and agreed with the literature. There was negligible impact on assay performance using FFPE tissues that had been archived for up to 30 years. Similarly, there was negligible impact of H&E staining, facilitating accurate visualization for scraping and assay of small focal areas of specific histology within a section.

Comparison of base-line and chemical-induced transcriptomic responses in HepaRG and RPTEC/TERT1 cells using TempO-Seq.

The utilisation of genome-wide transcriptomics has played a pivotal role in advancing the field of toxicology, allowing the mapping of transcriptional signatures to chemical exposures. These activities have uncovered several transcriptionally regulated pathways that can be utilised for assessing the perturbation impact of a chemical and also the identification of toxic mode of action. However, current transcriptomic platforms are not very amenable to high-throughput workflows due to, high cost, complexities in sample preparation and relatively complex bioinformatic analysis. Thus, transcriptomic investigations are usually limited in dose and time dimensions and are, therefore, not optimal for implementation in risk assessment workflows. In this study, we investigated a new cost-effective, transcriptomic assay, TempO-Seq, which alleviates the aforementioned limitations. This technique was evaluated in a 6-compound screen, utilising differentiated kidney (RPTEC/TERT1) and liver (HepaRG) cells and compared to non-transcriptomic label-free sensitive endpoints of chemical-induced disturbances, namely phase contrast morphology, xCELLigence and glycolysis. Non-proliferating cell monolayers were exposed to six sub-lethal concentrations of each compound for 24 h. The results show that utilising a 2839 gene panel, it is possible to discriminate basal tissue-specific signatures, generate dose-response relationships and to discriminate compound-specific and cell type-specific responses. This study also reiterates previous findings that chemical-induced transcriptomic alterations occur prior to cytotoxicity and that transcriptomics provides in depth mechanistic information of the effects of chemicals on cellular transcriptional responses. TempO-Seq is a robust transcriptomic platform that is well suited for in vitro toxicity experiments.

A trichostatin A expression signature identified by TempO-Seq targeted whole transcriptome profiling.

The use of gene expression signatures to classify compounds, identify efficacy or toxicity, and differentiate close analogs relies on the sensitivity of the method to identify modulated genes. We used a novel ligation-based targeted whole transcriptome expression profiling assay, TempO-Seq®, to determine whether previously unreported compound-responsive genes could be identified and incorporated into a broad but specific compound signature. TempO-Seq exhibits 99.6% specificity, single cell sensitivity, and excellent correlation with fold differences measured by RNA-Seq (R2 = 0.9) for 20,629 targets. Unlike many expression assays, TempO-Seq does not require RNA purification, cDNA synthesis, or capture of targeted RNA, and lacks a 3' end bias. To investigate the sensitivity of the TempO-Seq assay to identify significantly modulated compound-responsive genes, we derived whole transcriptome profiles from MCF-7 cells treated with the histone deacetylase inhibitor Trichostatin A (TSA) and identified more than 9,000 differentially expressed genes. The TSA profile for MCF-7 cells overlapped those for HL-60 and PC-3 cells in the Connectivity Map (cMAP) database, suggesting a common TSA-specific expression profile independent of baseline gene expression. A 43-gene cell-independent TSA signature was extracted from cMAP and confirmed in TempO-Seq MCF-7 data. Additional genes that were not previously reported to be TSA responsive in the cMAP database were also identified. TSA treatment of 5 cell types revealed 1,136 differentially expressed genes in common, including 785 genes not previously reported to be TSA responsive. We conclude that TSA induces a specific expression signature that is consistent across widely different cell types, that this signature contains genes not previously associated with TSA responses, and that TempO-Seq provides the sensitive differential expression detection needed to define such compound-specific, cell-independent, changes in expression.

A chemical-biological similarity-based grouping of complex substances as a prototype approach for evaluating chemical alternatives.

Comparative assessment of potential human health impacts is a critical step in evaluating both chemical alternatives and existing products on the market. Most alternatives assessments are conducted on a chemical-by-chemical basis and it is seldom acknowledged that humans are exposed to complex products, not individual substances. Indeed, substances of Unknown or Variable composition, Complex reaction products, and Biological materials (UVCBs) are ubiquitous in commerce yet they present a major challenge for registration and health assessments. Here, we present a comprehensive experimental and computational approach to categorize UVCBs according to global similarities in their bioactivity using a suite of in vitro models. We used petroleum substances, an important group of UVCBs which are grouped for regulatory approval and read-across primarily on physico-chemical properties and the manufacturing process, and only partially based on toxicity data, as a case study. We exposed induced pluripotent stem cell-derived cardiomyocytes and hepatocytes to DMSO-soluble extracts of 21 petroleum substances from five product groups. Concentration-response data from high-content imaging in cardiomyocytes and hepatocytes, as well as targeted high-throughput transcriptomic analysis of the hepatocytes, revealed distinct groups of petroleum substances. Data integration showed that bioactivity profiling affords clustering of petroleum substances in a manner similar to the manufacturing process-based categories. Moreover, we observed a high degree of correlation between bioactivity profiles and physico-chemical properties, as well as improved groupings when chemical and biological data were combined. Altogether, we demonstrate how novel in vitro screening approaches can be effectively utilized in combination with physico-chemical characteristics to group complex substances and enable read-across. This approach allows for rapid and scientifically-informed evaluation of health impacts of both existing substances and their chemical alternatives.

Feasibility of using microbeads with holographic barcodes to track DNA specimens in the clinical molecular laboratory.

We demonstrate the feasibility of using glass microbeads with a holographic barcode identifier to track DNA specimens in the molecular pathology laboratory. These beads can be added to peripheral blood specimens and are carried through automated DNA extraction protocols that use magnetic glass particles. We found that an adequate number of microbeads are consistently carried over during genomic DNA extraction to allow specimen identification, that the beads do not interfere with the performance of several different molecular assays, and that the beads and genomic DNA remain stable when stored together under regular storage conditions in the molecular pathology laboratory. The beads function as an internal, easily readable specimen barcode. This approach may be useful for identifying DNA specimens and reducing errors associated with molecular laboratory testing.

Medium- to high-throughput SNP genotyping using VeraCode microbeads.

Recent breakthroughs in multiplexed SNP (single nucleotide polymorphism) genotyping technology have enabled global mapping of the relationships between genetic variation and disease. Discoveries made by such whole-genome association studies often spur further interest in surveying more focused subsets of SNPs for validation or research purposes. Here we describe a new SNP genotyping platform that is flexible in assay content and multiplexing (up to 384 analytes), and can serve medium- to high-throughput applications. The Illumina BeadXpress platform supports the GoldenGate Genotyping Assay on digitally inscribed VeraCode microbeads to allow streamlined workflow, rapid detection, unparalleled data reproducibility and consistency. Thus, it is a highly valuable tool for biomarker research and validation, pharmaceutical development, as well as the development of molecular diagnostic tests.

Quantitative expression profiling of RNA from formalin-fixed, paraffin-embedded tissues using randomly assembled bead arrays.

Formalin-fixed, paraffin-embedded (FFPE) tissues represent an invaluable resource for gene expression analysis, as they are the most widely available materials for studies of human disease. However, degradation of RNA during tissue fixation and storage makes FFPE-derived RNAs hard to work with using conventional microarray technology. Most gene expression studies done using FFPE tissues rely on quantitative RT-PCR (qPCR), but this approach has had limited success because of the short cDNA templates available in these samples. In this chapter, we describe the DASL (cDNA-mediated annealing, selection, extension, and ligation) Assay, a flexible, sensitive, and reproducible gene expression profiling system for parallel analysis of several hundred mRNA transcripts on the BeadArray platform. This technology is especially useful for determining cancer prognosis or therapy response, because it allows not only prospective studies but also retrospective analyses. Using the DASL Assay, gene expression analyses can be performed on routinely stored tumor specimens from patients with known outcomes.

Spatial analysis of arabidopsis thaliana gene expression in response to Turnip mosaic virus infection.

Virus-infected leaf tissues comprise a heterogeneous mixture of cells at different stages of infection. The spatial and temporal relationships between sites of virus accumulation and the accompanying host responses, such as altered host gene expression, are not well defined. To address this issue, we utilized Turnip mosaic virus (TuMV) tagged with the green fluorescent protein to guide the dissection of infection foci into four distinct zones. The abundance of Arabidopsis thaliana mRNA transcripts in each of the four zones then was assayed using the Arabidopsis ATH1 GeneChip oligonucleotide microarray (Affymetrix). mRNA transcripts with significantly altered expression profiles were determined across gradients of virus accumulation spanning groups of cells in and around foci at different stages of infection. The extent to which TuMV-responsive genes were up- or downregulated primarily correlated with the amount of virus accumulation regardless of gene function. The spatial analysis also allowed new suites of coordinately regulated genes to be identified that are associated with chloroplast functions (decreased), sulfate assimilation (decreased), cell wall extensibility (decreased), and protein synthesis and turnover (induced). The functions of these downregulated genes are consistent with viral symptoms, such as chlorosis and stunted growth, providing new insight into mechanisms of pathogenesis.

A regulator of Dscam mutually exclusive splicing fidelity.

The Down syndrome cell adhesion molecule (Dscam) gene has essential roles in neural wiring and pathogen recognition in Drosophila melanogaster. Dscam encodes 38,016 distinct isoforms via extensive alternative splicing. The 95 alternative exons in Dscam are organized into clusters that are spliced in a mutually exclusive manner. The exon 6 cluster contains 48 variable exons and uses a complex system of competing RNA structures to ensure that only one variable exon is included. Here we show that the heterogeneous nuclear ribonucleoprotein hrp36 acts specifically within, and throughout, the exon 6 cluster to prevent the inclusion of multiple exons. Moreover, hrp36 prevents serine/arginine-rich proteins from promoting the ectopic inclusion of multiple exon 6 variants. Thus, the fidelity of mutually exclusive splicing in the exon 6 cluster is governed by an intricate combination of alternative RNA structures and a globally acting splicing repressor.

Illumina universal bead arrays.

This chapter describes an accurate, scalable, and flexible microarray technology. It includes a miniaturized array platform where each individual feature is quality controlled and a versatile assay that can be adapted for various genetic analyses, such as single nucleotide polymorphism genotyping, DNA methylation detection, and gene expression profiling. This chapter describes the concept of the BeadArray technology, two different Array of Arrays formats, the assay scheme and protocol, the performance of the system, and its use in large-scale genetic, epigenetic, and expression studies.

Two-dimensional transcriptome profiling: identification of messenger RNA isoform signatures in prostate cancer from archived paraffin-embedded cancer specimens.

The expression of specific mRNA isoforms may uniquely reflect the biological state of a cell because it reflects the integrated outcome of both transcriptional and posttranscriptional regulation. In this study, we constructed a splicing array to examine approximately 1,500 mRNA isoforms from a panel of genes previously implicated in prostate cancer and identified a large number of cell type-specific mRNA isoforms. We also developed a novel "two-dimensional" profiling strategy to simultaneously quantify changes in splicing and transcript abundance; the results revealed extensive covariation between transcription and splicing in prostate cancer cells. Taking advantage of the ability of our technology to analyze RNA from formalin-fixed, paraffin-embedded tissues, we derived a specific set of mRNA isoform biomarkers for prostate cancer using independent panels of tissue samples for feature selection and cross-analysis. A number of cancer-specific splicing switch events were further validated by laser capture microdissection. Quantitative changes in transcription/RNA stability and qualitative differences in splicing ratio may thus be combined to characterize tumorigenic programs and signature mRNA isoforms may serve as unique biomarkers for tumor diagnosis and prognosis.

Tobamovirus infection is independent of HSP101 mRNA induction and protein expression.

Heat shock protein 101 (HSP101) has been implicated in tobamovirus infections by virtue of its ability to enhance translation of mRNAs possessing the 5'Omega-leader of Tobacco mosaic virus (TMV). Enhanced translation is mediated by HSP101 binding to a CAA-repeat motif in TMV Omega leader. CAA repeat sequences are present in the 5' leaders of other tobamoviruses including Oilseed rape mosaic virus (ORMV), which infects Arabidopsis thaliana. HSP101 is one of eight HSP100 gene family members encoded by the A. thaliana genome, and of these, HSP101 and HSP98.7 are predicted to encode proteins localized to the cytoplasm where they could potentially interact with TMV RNA. Analysis of the expression of the HSP100s showed that only HSP101 mRNA transcripts were induced significantly by ORMV in A. thaliana. The induction of HSP101 mRNA was also correlated with an increase in its protein levels and was independent of defense-related signaling pathways involving salicylic acid, jasmonic acid, or ethylene. A. thaliana mutants lacking HSP101, HSP98.7, or both supported wild-type levels of ORMV replication and movement. Similar results were obtained for TMV infection in Nicotiana benthamiana plants silenced for HSP101, demonstrating that HSP101 is not necessary for efficient tobamovirus infection.

Salicylic acid-dependent expression of host genes in compatible Arabidopsis-virus interactions.

Plant viruses elicit the expression of common sets of genes in susceptible hosts. Studies in Arabidopsis (Arabidopsis thaliana) and tomato (Lycopersicon esculentum) indicate that at least one-third of the genes induced in common by viruses have been previously associated with plant defense and stress responses. The genetic and molecular requirements for the induction of these stress and defense-related genes during compatible host-virus interactions were investigated with a panel of Arabidopsis mutant and transgenic plants defective in one or more defense signaling pathways. pad4, eds5, NahG, npr1, jar1, ein2, sid2, eds1, and wild-type Columbia-0 and Wassilewskija-2 plants were infected with two different viruses, cucumber mosaic virus and oilseed rape mosaic virus. Gene expression was assayed by a high-throughput fiber-optic bead array consisting of 388 genes and by RNA gel blots. These analyses demonstrated that, in compatible host-virus interactions, the expression of the majority of defense-related genes is induced by a salicylic acid-dependent, NPR1-independent signaling pathway with a few notable exceptions that did require NPR1. Interestingly, none of the mutant or transgenic plants showed enhanced susceptibility to either cucumber mosaic virus or oilseed rape mosaic virus based on both symptoms and virus accumulation. This observation is in contrast to the enhanced disease susceptibility phenotypes that these mutations or transgenes confer to some bacterial and fungal pathogens. These experimental results suggest that expression of many defense-related genes in compatible host plants might share components of signaling pathways involved in incompatible host-pathogen interactions, but their increased expression has no negative effect on viral infection.

Quantitative gene expression profiling in formalin-fixed, paraffin-embedded tissues using universal bead arrays.

We recently developed a sensitive and flexible gene expression profiling system that is not dependent on an intact poly-A tail and showed that it could be used to analyze degraded RNA samples. We hypothesized that the DASL (cDNA-mediated annealing, selection, extension and ligation) assay might be suitable for the analysis of formalin-fixed, paraffin-embedded tissues, an important source of archival tissue material. We now show that, using the DASL assay system, highly reproducible tissue- and cancer-specific gene expression profiles can be obtained with as little as 50 ng of total RNA isolated from formalin-fixed tissues that had been stored from 1 to over 10 years. Further, tissue- and cancer-specific markers derived from previous genome-wide expression profiling studies of fresh-frozen samples were validated in the formalin-fixed samples. The DASL assay system should prove useful for high-throughput expression profiling of archived clinical samples.

A versatile assay for high-throughput gene expression profiling on universal array matrices.

We report a flexible, sensitive, and quantitative gene-expression profiling system for assaying more than 400 genes, with three probes per gene, for 96 samples in parallel. The cDNA-mediated annealing, selection, extension and ligation (DASL) assay targets specific transcripts, using oligonucleotides containing unique address sequences that can hybridize to universal arrays. Cell-specific gene expression profiles were obtained using this assay for hormone-treated cell lines and laser-capture microdissected cancer tissues. Gene expression profiles derived from this assay were consistent with those determined by qRT-PCR. The DASL assay has been automated for use with a bead-based 96-array matrix system. The combined high-throughput assay and readout system is accurate and efficient, and can cost-effectively profile the expression of hundreds of genes in thousands of samples.

Dilated cardiomyopathy caused by tissue-specific ablation of SC35 in the heart.

Many genetic diseases are caused by mutations in cis-acting splicing signals, but few are triggered by defective trans-acting splicing factors. Here we report that tissue-specific ablation of the splicing factor SC35 in the heart causes dilated cardiomyopathy (DCM). Although SC35 was deleted early in cardiogenesis by using the MLC-2v-Cre transgenic mouse, heart development appeared largely unaffected, with the DCM phenotype developing 3-5 weeks after birth and the mutant animals having a normal life span. This nonlethal phenotype allowed the identification of downregulated genes by microarray, one of which was the cardiac-specific ryanodine receptor 2. We showed that downregulation of this critical Ca2+ release channel preceded disease symptoms and that the mutant cardiomyocytes exhibited frequency-dependent excitation-contraction coupling defects. The implication of SC35 in heart disease agrees with a recently documented link of SC35 expression to heart failure and interference of splicing regulation during infection by myocarditis-causing viruses. These studies raise a new paradigm for the etiology of certain human heart diseases of genetic or environmental origin that may be triggered by dysfunction in RNA processing.

Expression of an active tobacco mitogen-activated protein kinase kinase kinase enhances freezing tolerance in transgenic maize.

Cold acclimation is the major process that prepares plants for freezing tolerance. In addition to extensive transcription regulation by cold-inducible master transcription factors, oxidative stress signaling has been postulated to play a role in freezing tolerance. Activation of oxidative signaling through the expression of an active mitogen-activated protein kinase kinase kinase provided benefits in transgenic tobacco at freezing temperature bypassing cold acclimation. Because involvement of the mitogen-activated protein kinase cascade in oxidative stress signaling is evolutionarily conserved in eukaryotes from yeast to mammals, we tested the effect of expressing a heterologous tobacco mitogen-activated protein kinase kinase kinase (Nicotiana PK1), which can mimic H(2)O(2) signaling, in a major cereal crop. We demonstrate that low-level but constitutive expression of the Nicotiana PK1 gene enhances freezing tolerance in transgenic maize plants that are normally frost sensitive. Our results suggest that a new molecular approach can be designed to genetically enhance freezing tolerance in important crops.

Profiling alternative splicing on fiber-optic arrays.

The human transcriptome is marked by extensive alternative mRNA splicing and the expression of many closely related genes, which may be difficult to distinguish using standard microarray techniques. Here we describe a sensitive and specific assay for parallel analysis of mRNA isoforms on a fiber-optic microarray platform. The method permits analysis of mRNA transcripts without prior RNA purification or cDNA synthesis. Using an endogenously expressed viral transcript as a model, we demonstrated that the assay readily detects mRNA isoforms from as little as 10-100 pg of total cellular RNA or directly from a few cells. Multiplexed analysis of human cancer cell lines revealed differences in mRNA splicing and suggested a potential autocrine mechanism in the development of choriocarcinomas. Our approach may be useful in the large-scale analysis of the role of alternative splicing in development and disease.